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Abstract 

It is the first time that the study of three-body baryonic B decays offers an independent deter- 
mination of the nucleon form factors for the time-like four momentum transfer (t > 0) region, such 
as G P j^{t) of vector currents and those of axial ones. Explicitly, from the data of B° — ► npD* + 
and A.pTi + we find a constant ratio of G n M (t) /G p M (t) = -1.3 ± 0.4, which supports the FENICE 
experimental result. The vector and axial-vector form factors of pn, pp and nn pairs due to weak 
currents are also presented, which can be tested in future experiments. 

PACS numbers: 13.40.Gp, 28.20.-v, 13.25.Hw, 14.20.Dh 
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The nucleon form factors still attract attentions as they reveal the hadron structures and 
play important roles in any scattering or decaying processes involving baryons jjj. These 
form factors depend on the four momentum transfer (t), which is either space-like (t < 0) or 
time-like (t > 0). The behaviors of the form factors versus t have been extensively studied 



m various 



QCD models 
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131 ] . such as perturbative QCD 



(PQCD), chiral perturbation theory (ChPT), vector meson dominance (VMD) approach and 
dispersion relation (DR) method. Experimentally, accurate electromagnetic data for vector 
form factors with the space-like momentum transfer have become abundant [jjj], whereas 
the data on axial- vector form factors are available only for the space-like region with \t\ < 1 



GeV 2 in the neutrino-nucleon scattering [l|. In particular, the time-like electromagnetic form 

factors of the proton have been extracted from e + e~ — > pp (pp — > e + e~) 15|, but only a 

:'ew data points have been collected for those of the neutron by the FENICE Collaboration 
~l 

161 ] . Currently, due to experimental difficulties, there are no data on the time-like axial 
structures, induced from the weak currents due to W and Z bosons. Moreover, t 
some inconsistency between the measurement and theory {2], 3], 0, Q, E3, Q, [§], [9, 10, 
for the nn data, unlike the pp case, which seems to be well understood by the theoretical 
calculations. Clearly, more theoretical studies as well as precise experimental measurements 
on the nucleon form factors are needed to improve our understanding of strong interactions. 

In this paper, we shall show that the three-body baryonic B decays of B — > BB'M, such 
as B° — > npD* + and B° — > Apir + , can provide valuable information on the nucleon form 
factors. In general, the three-body baryonic B decays involve time-like form factors from 
vector, axial vector, scalar and pseudoscalar currents, respectively. In the scale of rrib ~ 4 
GeV, the PQCD is suitable for us to systematically examine not only the form factors of 



rere exists 
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vector currents but also those o 
series of works in Refs. 
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axial vector ones. We note that the PQCD approach in a 
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231 ] has been developed reliable tool to 



explain the experimental data on the baryonic B decays. 

In the wildly used factorization method 2J,l25[, which splits the four quark operators into 
two currents by the vacuum insertion, there are three types of three-body baryonic B decays: 
Type I is for the decay in which a meson is transformed from B together with an emitted 
baryon pair; Type II is for the mode in which a baryon pair is transited from B together 
with an ejected meson; and Type III is the mixture of Types I and II. With the factorization 
method, the decay amplitudes for Types I and II are proportional to (BB'| J^|0) (M \ J%\B) 



2 



and {M\Jl 



|0)(BB'| J 2 \B), respectively. For the present measured modes, for instance, B 



28| and B~ 



npD* + M and B° -> Apir + |27| belong to Type I, while B° -> ppD^° 
KpJ/^ 29] are classified as Type II, whereas B~ — > ppK^~ , 5° — > ppK$, B~ — > ppix~ 
and -B - — > AAA' - are of Type III. Although the decay modes of Types I and III are 
of our current interest, those of Type III are inevitably affected by the uncertainties of the 
B — > BB' transition form factors. In this study, we shall concentrate on the Type I modes 
of B° -> npD* + and B° -> Apn + . 



by 
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With the effective Hamiltonians 32| at the quark level, the decay amplitudes are given 
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22] 



G f 



V cb V: dai (np\(du) v . A \0)(D* + \(cb) v _ A \B°) 



(1) 



,4(S ^Ap7r+) = -||(K fe K>i- W;a4)(Ap|(^) y _ A |0)(7r + |(^) y _ A | J B°) 
+V tb V; s 2a e (Ap\(su) s+P \0}(n + \(ub) s ^ P \B } 



(2) 



where is the Fermi constant, V qiqj are the CKM matrix elements, (qiqj)v-A = 9i7/i(l 
75) 9j, (QiQj)s±p = ± 75)^ and a* (z = 1, 4, 6) are given by 



ax 



eff 1 



CI4 



+ 



A 7 , 



a 6 



+ 



(3) 



with c 



= 1, 2, ■ ■ • , 6) being the effective Wilson coefficients (WCs) shown in Refs. 
321 ] and N**' the effective color number. Here, we have used the generalized factorization 
method with the non-factorizable effect absorbed in A 7 "^ ■ For the matrix elements of 
(ir + \(ub) v _ A \B°) and (D* + \(cb) v _ A \B°), we use the results in Refs. 
like baryonic form factors, we have 
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3311 - For the time- 



(BB'|g i7 ^-|0) = u(p B )\F 1 (t) llM + 



u(p B )UF 1 (t) + F 2 (t)}^ + 



m B + mg 



(Pb ~Pb)h }v(Pb 



J B) , 



(BB / |g i 7 M 7 5 qf J -|0) = u(p B ) i #a(*)7m + 



h A (t) 



?riB + mi 



■(Pb>+Pb)h >l5v{p 



J B) , 



(BB'l^-10) = fs(t)u(p B )v(p B ,) , 
(BB'|ft 75 ty-|0) = gp(t)u(pB)75v(p B ,) , 



(4) 



where the four momentum transfer in the time-like region is t = (p B + P B <) 2 , Qi = u , d an d 
s, and Fi, F 2 , g A , h A , f s and gp are the form factors. 



In this paper, we will study the form factors in Eq. (j4j) based on the experimental data 
in the baryonic B decays. We begin by defining the baryonic form factors of F\ and qa by 



(B|JT|B / >=fi(Ri 



*i(*)7/i + &i(<)7M7B 



v(pi 



(5) 



where Jf™ 



QgQlfi Q an d t = (pb — Pb 1 ) 2 - Note that F 2 and h A are not included in Eq. (JS]) 
due to the helicity conservation. To exhibit the chirality or helicity, we rewrite Eq. §5§ as 



(B T+ ;|J^|BVi>=fi(PB) 



7m" 



75 



1 



7m" 



75-1 



G (t) 



u(pi 



(6) 



2 v ' "* 2 

where |B-|- + j_) = |B|) + |Bj_) respects both flavor SU(3) and spin SU(2) symmetries, e.g., 
| Pl) = ^/l/18[u|U|<i| + iijU|d| — 2u|«|rf| + permutations]. In Eq. (jSJ), G T (t) and G^(t) 
represent the right-handed and left-handed form factors, which can be further decomposed 

as 



(7) 



id) k \ ) 

where the constants e , and ei are defined by 



J(4) 



( B ta)IQiil B m)) ' 4 a) = ( B Ta)IQnl B Ta)) 



l|l"T(i)/ ' 



(8) 



respectively, with Q 



(II) 



In Eq. (jSJ), the summation is over the charges carried 



by the valence quarks (i — 1,2, 3) in the baryon with helicities parallel (||) and anti-parallel 
(||) to the baryon spin directions of (j, j). Since are the form factors accompanied 

by the (anti-)parallel hard-scattering amplitudes with baryon wave functions, based on the 



QCD counting rules in the PQCD 
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351 ] . they can be expressed by 



G\\(t) 



t 2 



"4> 



1 -7 



Gj>(t) 



Ctt 



f 1 



1 



(9) 



where 7 = 2.148, A = 300 MeV 



6[ and Cm tt are parameters to be determined. We remark 
that the asymptotic formulas in Eq. (Q are exact only when t is large. For smaller t, such 
as when t being close to the two-nucleon threshold, higher power corrections are expected 



36( | . The corrections will be averaged into the errors of C,, tt in our data fitting. From Eqs. 



©- (ED, we have 



F 1 (t) = (el + e^G^ + iel + e^G-fi) 



g A (t) = (eJ- e J)G||(0 + (4-4)CT|(t). 



(10) 



te region, the time-like form factors 



Although the above equations are derived in the space-lil 

can be easily written via the crossing symmetry {(], [37], which transforms the particle in 
the initial state to its anti-particle in the final state and reverses its helicity. However, in 
general, the values of the time-like G\\(t) and G-^(t) form factors are complex numbers unlike 
the space-like ones for which there is a time reversal symmetry between initial and final 
states. 



In Eq. (|4]), F 2 is suppressed by l/^tln^/A 2 ,]) in comparison with F 1 [38|, |39| and therefore 
can be safely ignored, while g P is found to be related to fs as j2o| 

9p = fs , (11) 
and fs and Ha are deduced from equation of motion, given by 

M( ) = J m " + t m ^ gA{t) . (12) 

Thus, once we figure out F\ and gA, all other form factors fs, Ha and gp will be determined 
in terms of Eqs. (TTTT) and (|T2|) . For the electromagnetic current J^ m = ^u^^u — ^d'j^d, from 
Eq. ffTOl) it is clear that = since e^-y = ej - . Furthermore, from Eqs. dl])-® we have 

G p M (t) = Ff{t)+Ff{t) ~ FfJ m) (t) = G||(t), 

g^w = Fr(t)+Frw * ^ n ( L)(t) = _^ + 3j^, (is) 



where we have neglected the small F 2 terms comparing with those of F^ [38|, l39 |. 



Similarly, for the Z coupled current = \u r y ll {— ^ i )w — \d^^{^— —)d — sin 2 9wJ^ m , we use 

tZ 



the same form factors defined in Eqs. (0) and ([6]) by replacing J^ m with and we get 
= 2 ~ 3 f 2 ^ Gn(t) - -^(t) , <^(t) = |C||(t) + , 



6 ||W , 3 6 

= - r — G\\(t) + 1 + 2 f 2 ^ %(t), ^(t) = ~G\\{t) - iGtft) .(14) 



-2-sin 2 #w ,. N l + 2sin 2 flw „ a _ 2 _ s 1 



3 liw 6 ||v " / 3 MW g 

Since the behaviors of G\\(t) and Gjr(i) have been given in Eq. (jUj), what we shall do next is 
to fix the parameters C\\ and Cry in terms of 

^ Ap "w = yfcn(t) , ^(t) = J|g,i(*) , (is) 



derived from Eqs. ([I])-® with the data in B° — > npD* + and B° — > Ap7i + . 

Before performing the numerical analysis, we would like briefly discuss the generalized 
factorization method. It is known that the factorization method 
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251 ] suffers from several 



possible hadron uncertainties, such as those from the nonfactorizable effect, annihilation 
contribution and final state interaction. To describe these uncertainties, we take the decay 
of B° — > npD* + as an example, while those for B° — > Ap7T + can be treated in a similar 
manner. The amplitude of B° — > npD* + from the color suppressed operator is given by 



c 2 ft (npD* + \(d a U{3) V -A(cpb a ) V -A\B ) 



-,eff 



(np\(du) v _ A \0)(D* + \(cb) v _ A \B o ) + -*-(npD* + \(d\ a u) v _ A (c\ a b) v _ A \B°) ,(16) 



where 5ppi5 aa i = 5f3 a 5 a >f3>/N c + A^A^/2 has been used to deal with color index a (/?) and 
the second term on the right-hand side is the so-called nonfactorizable effect. Although 
the nonfactorizable effect cannot be directly and unambiguously determined by theoretical 
calculations, in the generalized factorization method 3^] , this contribution can be absorbed 
in the effective color number N*** running from 2 to oo in Eq. Q. The amplitude of the 
annihilation contribution is given by 

Aan(B° ^npD* + ) = ^V cb V: d a 2 (npD* + \(cu) v . A \0}(0\(db) v . A \B } , (17) 

V2 

where a 2 is around (0.1 — 0.01)ai, which is suppressed. In addition, based on the power ex- 
pansion of 1/t in the PQCD approach, A an (B° npD* + ) oc 1/t 3 , which is much suppressed 
than A(B° -> npD* + ) as t ~ ml [40J . For the final state interaction, the most possible source 
is via the two-particle rescattering to the baryon pair, such as B° — > MiM 2 D* + — > npD* + 
with Mi 2 representing meson states. However, such processes would shape the curve asso- 



ciated with the phase spaces in the decay rate distributions 41 
excluded in the charmless baryonic B decay experiments 
In our numerical analysis, we take 



42 



43], which have been 



1.17, -0.37, 0.0246, -0.0523, 0.0154-0.066) 



32 



effr ■ 



44 



l,2,---,6) 
45j, m u (m b ) =3.2 MeV, mJm h ) = 90 



MeV |32| and m h {m b ) = 4.2 GeV j46j, and the weak phase 7 = 59.8° ± 4.9° [47j. For the 
experimental data in the B decays, we use |26l. l27| 



Br(B° -> npD*^ 
Br(B° -> Ap7r A 



(14.5 ±4.3) x 10" 4 , 
(3.29 ±0.47) x 10~ 6 



(18) 



and the other available data in Ref. [27j, such as spectrum vs. invariant mass and angular 
distributions. 

Based on the x 2 fitting, we obtain 



Y 



'■/dof = 0.7 , l/N e J f = 0.2 ± 0.2, 



(19) 



where dof denotes the degree of freedom. It is clear that x 2 /dof = 0.7 presents a reliable fit, 
while A^f ~ 2.5 — oo means a limited nonfactorizable effects with the small contributions 
from the annihilation and final state interaction. We stress that if the hadronic uncertainties 
were large, N^* from 2 to oo would not be account for the data. The coefficients of C\\ and 
Cjj in Eq. IQ are found to be 



which lead to 



Cm = 83.7 ±5.7 Gel/ 4 and C Tl 



-246.3 ±92.1 Gel/ 4 



(20) 



%(t)/G,|(t) 



-2.9 ± 1.1 



(21) 
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since their imaginary parts are expected 
371 ] as well as the result in the DR method 



Here, we have assumed that both C\\ and Cry are rea 
;o be small based on the argument in Refs. 

As seen in Fig. [1^, the fitted values of G 7 ^ (t) extracted from the B decays are in agreement 
with the FENICE data fl6|] by the assumptions of |G^-| = \G%\ and \G%\ = 0. We note that 
G~%{i) = F™ n (t) ± jrrt-F™ is the neutron electric form factor. Clearly, in our calculation 
based on the QCD counting rule, G% ~ G%j. 

From Eqs. (191. ([TH]) and (|2Q|1 , we get 



G n M (t)/G p M (t) = -1.3 ±0.4 



(22) 



which supports the measurements in Refs. 15| and 16j. We note that in Eq. (122]) the ratio 
is a constant due to the same power expansions in Eq. (Q and the minus sign is necessary 
in order to match the measured baryonic B decay branching ratios, which also confirms the 



theoretical result in Refs 



34 



351 ] . Moreover, the puzzle of IG^^t) / G p M (t)\ > 1 is also 



solved as indicated in Eq. (1221) . We note that the ratio G^(t) / G p M {t) is predicted to be —2/3 
and —1/2 in the QCD counting Q] and sum rules {3J, respectively, while the DR 4, 5, 6, 7| 
and VMD 0, methods yield only half of the values indicated by the data points. It is 




PIG. 1: Fo rm factors of (a, <*(«), (o) and (c> ttW with toe Hm e* four _ 
transfer i, where the star and triangle symbols represent the FENICE data [lfj with the assump- 
tions of = \G%\ and |G^| = 0, and the solid, dash and dotted curve stand for (b) Ff n (t), 
Ff {z) {t) and F™ z) {t) and (c) gf(t), gf (z) (t) and gf {z) (t), respectively. 

interesting to see that if G\\(t) / Gjr(t) = —1 instead of the fitted one in Eq. (l2Tj) . from Eq. 



Therefore, our result on the ratio in 



(USD G^{t)/G p M {t) = -2/3 is recovered as in Ref. 
Eq. (I2T1) could help us to improve the PQCD calculations. 

Since we have related all the nucleon form factors as in Eqs. ( [TBI . ( |T4|) and ( |T5l) . we find 

Ff{t) = (193.7 ± 31.6)G p (t) , gf(t) = (29.5 ± 31.6)G p (t) , 
Fff z) (t) = (78.4 ± 15.6)G P (0 , ^ z) (t) = (14.8 ± 15.8)G p (t) , 

= (-121.1 ± 22.5)G p (t) , gf {z) (t) = (-14.8 ± 15.8)G p (t) , (23) 

where G p (t) = l/t 2 [ln(t/Ag)]" 7 . The central values of F^' and g^' in Eq. ([23]) are shown 
in Figs. [T]d and [It, respectively. The valid ranges for these time-like form factors are the 
same as those in the three-body baryonic B decays of B — > BB'M, i.e., (hib + rn&>) 2 — 
4 GeV 2 < t < [m B - m M ) 2 ~ 16 - 25 GeV 2 . 

In sum, we have shown that the study of the measured three-body baryonic B decays 
of B° -> npD* + and B° Apn + leads to G n P {t) /G p t {t) = -1.3 ± 0.4, which supports the 
FENICE measurement. The minus sign for the ratio Gp(t)/Gf(t), given by the previous 
theoretical calculations, has been enforced to fit the B decay data. We have pointed out 
that our fitted value for the ratio of G\\(t) and Grr(t) may be useful for us to perform various 
QCD calculations on e + e~ — > nn. We have also predicted the time-like vector and axial 
vector nucleon form factors induced from the weak currents, such as Ff N '(t) and g^ N ' (t) 
(N, N' =p and n). 

Finally, we remark that apart from the use of B° — > pnD* + and B° — > ApTT + , there are 
more decays directly connecting to the time-like form factors, such as B° — > pn(D + , p + , ir + ) 

8 



and B° — > App + as well as the corresponding charged B modes, which are within the 
accessibility of the current B factories at KEK and SLAC. It is clear that as more and 
more data available from current and future B factories, the nucleon form factors can 
be further constrained and determined. Moreover, the new measurements in e + e~ — > nn 
are progress in DA«E at FVascati fl and p.annin g in PANDA and PAX at GSI 
[3]. As for the weak nucleon form factors, since the scattering of e + e~ at BABAR is at 
the mB scale, Ff N '(t) and g^ N ' (t) (N,N' = p and n) can be studied via the left-right 



helicity asymmetry [36| of A PV = {daR—da L )/{da R +daL) as in the SAMPLE experiment [l|. 
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